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Abstract 

We study Wilson loops at non-zero temperature and extract the static quark potential from them. 
The extracted potentials are larger than the singlet free energies and do not show screening for 
T < 190 Me V. 



1. Introduction 

Quarkonium suppression was proposed by Matsui and Satz as a signature of formation of 
deconfined medium in heavy ion collisions uJ]. The basic idea behind this proposal is that color 
screening in the deconfined medium will modify the heavy quark potential, eventually leading 
to the dissolution of the heavy quarkonium states. The problem of dissolution of quarkonium 
\ states at high temperatures could be formulated in terms of spectral functions. Early attempts to 

• calculate spectral functions on the lattice have been presented in Refs. However, extraction 

of the spectral functions from lattice results on Euclidean correlation functions is quite difficult 
Jl] and one should also be careful with cutoff effects in the spectral functions extracted from the 
lattice 101 ■ Furthermore, the Euclidean time correlators may not be sensitive to the melting of 
the bound states at high temperatures due to the fact that the Euclidean time extent is limited 
to < l/(2T) (see e.g. discussions in Ref. jsj]). The effective field theory framework for heavy 
quark bound states, namely pQNRCD could be a useful tool for calculating quarkonium spectral 
functions [6J. The effective field theory approach allows to rigorously define the concept of 
the static quark anti-quark potential both at zero and non-zero temperatures. One of the main 
outcomes of the effective field theory analysis is the finding that at non-zero temperature the 
potential has also an imaginary part, which has important consequences for the dissolution of 
the quarkonium states. While pNRQCD is formulated in the weak coupling framework it is 
possible to extend it to the non-perturbative regime. For example, if the binding energy is the 
smallest scale in the problem all the other scales, like the thermal scales, the inverse size of the 
bound state and Aqcd can be integrated out. In this case the potential should be determined 
non-perturbatively and is identical to the energy of a static QQ pair If one further neglects 
the dipole interactions one gets the generalization of the simple potential model to the case of 
high temperatures However, one still needs to specify the potential. In the past model 
considerations based on lattice calculations of the so-called singlet free energy have been used 
(see e.g. discussion in Ref. liTJ). In Ref. |i8J it has been suggested to extract the energy of a static 
QQ pair using the spectral decomposition of the temporal Wilson loops at non-zero temperature. 
In this contribution we attempt to extract the static quark anti-quark energy in 2-1-1 flavor QCD 
based on this idea. 
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2. Numerical results 



In lattice QCD calculations the static QQ energy is extracted from Wilson loops Wir, t). 
At large Euclidean time separations the exponential decay of the Wilson loops is governed by 
the static energy or potential, W(r, t) ~ exp(-y(r)T). More generally one can write a spectral 
decomposition for the Wilson loops |18|] 

du)a-(r,u))e-'^\ (1) 

oo 

At zero temperature the spectral function is proportional to d(ix) - V(r)) plus a sum of delta 
functions corresponding to the excited states (hybrid potentials). At non-zero temperature the 
delta function becomes a Lorentzian with the width related to the imaginary part of the potential. 
In Ref. 1 8] the potential at non-zero temperature was extracted by inverting Eq. ([T]i via the 
maximum entropy method (MEM) which gives reasonably accurate determination of the real 
part of the potential. At the same time it is hard to estimate the imaginary part of the potential. 
The problem is similar to the problem of extracting meson spectral functions |2, 3], where the 
width of the bound state peaks is mostly an artifact of MEM. Another problem that appears in 
the calculation of the potential is that the Wilson loops become noisy at large spatial separations 
r. To deal with this problem smeared gauge fields are used in spatial links when constructing 
Wilson loops on the lattice. Alternatively, one can fix the Coulomb gauge and calculate the 
correlation functions of two temporal Wilson lines separated by distance r without connecting 
them by spatial links |9]. At zero temperature where one only interested in the energy levels both 
choices are equally good, and merely correspond to different choices of static meson interpolating 
operators. The same should be true at non-zero temperature provided the imaginary part of 
the potential is not too large. We have calculated the correlation functions of temporal Wilson 
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Figure 1 : The correlation function of Wilson lines as function of the time extent t calculated for 48''' X 12, yS = 7.5 (left) 
and 24^ X 6,/3 = 6.8 (right). 

lines in 2-1-1 flavor QCD using the Highly Improved Staggered Quark (HISQ) action fld\ and 
improved gauge action with the physical value of the strange quark mass and light quark masses 
corresponding to pion mass of 160 MeV in the continuum limit. The calculations have been 
performed on 48-^ x 16 and 48-^ x 12 lattices for bare gauge coupling /3 - lO/g^ = 7.5 as well 
as on 24^ x 6 lattices at various gauge cou plin gs that have been used in the study of the chiral 
and deconfinement transition by HotQCD jl ill . As in Ref. ifllll the lattice spacing was set by 
the ri scale extracted from the static potential and using the value ri - 0.3106 fm. We also 
used the renormalization constants for the Wilson line correlators determined in Ref. Ill 111 . It 
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should be noted that for tT - I the correlator gives the so-called singlet free energy (H]. The 
temperature dependence of the singlet free energy obtained in our calculations is very similar to 
the temperature dependence of the singlet free energy obtained earlier with the p4 action lfl3i[l4l] . 
In Fig. [1] we show our numerical results for the correlator of Wilson lines as function of the 
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Figure 2: The static quark anti-quark potential extracted from 24' X 6 lattices at different temperature compared to the 
zero temperature result as well as to the singlet free energy. 

Euclidean time extent r for different distances calculated on 48-^ x 12 lattice at/? = 7.5 and 24^ x6 
lattice at yS = 6.8. These gauge couplings correspond to temperatures 300 MeV and 320 MeV, 
respectively. As one can see from Fig. [T]the r-dependence of the correlators is consistent with 
simple exponential decay, except close to ^ 1, where the correlators show a slight increase. 
This increase is due to the contribution of a backward propagating state ~ exp(-£'/,flei(l/r - r)), 
that arises from the fact that static quarks propagate in gauge field background that is periodic in 
T. Similar behavior has been observed in the Wilson loop calculations at non-zero temperature 
in pure gauge theory [8J as well as in full QCD calculations of bottomonium spectral functions 
within the non-relativistic formulation OlSll . For the largest two values of A^^ considered here the 
backward propagating state does not cause any problem and we get stable results for the potential 
by performing single exponential fits in the r-interval around the mid-point tT - 1/2. However, 
the results are quite noisy for rT" > 1 in this case. Thus to explore the potential at larger distances 
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we use 24^ x 6 lattices for which statistical errors are small. Here the results are sensitive to how 
the fits are done. We performed three type of fits. First we used only tT - 1/3 and tT - 1/2 to 
extract the potential. Then we determined the backward propagating contribution by performing 
fits with tT - \ and tT - Ijli. Subtracting the backward propagating contribution from the 
correlator we extracted the potential using tT - 1/2 and tT - Ijli points. This gives our central 
value of the extracted potential. Finally, we performed single exponential fits for tT - Ijl and 
tT = 2/3 which gave us the lower value of the potential. 

Our numerical results for the static quark anti-quark potential at different temperatures are 
shown in Fig. |2l The errors for the temperature dependent potential shown in the figure are 
mostly systematics and are estimated as described above. We also compare the static quark 
anti-quark potential with the corresponding zero temperature results as well as with the singlet 
free energy. For the lowest temperature both the potential and the singlet free energy agree with 
the zero temperature result. For T - 178MeV the singlet free energy is very different from 
the zero temperature potential, while the difference is small for the finite temperature potential. 
Furthermore, the in-medium potential does not show screening at this temperature. Screening 
effects become apparent in the potential at T = 194MeV and happen at distances of about Ifm. 
At the same temperature the screening effects in the singlet free energy set in at distance of about 
0.6fm. At higher temperatures screening effects in the potential set in at smaller and smaller 
distances and the difference between the potential and the singlet free energy becomes smaller. 
This is ejcpected as at very high temperatures the singlet free energy should be equal to the 
potential At all temperatures the potential is larger than the singlet free energy, i.e. it seems 
to approach the singlet free energy from above. 

3. Conclusion 

We have calculated the correlation functions of Wilson lines at non-zero temperature for dif- 
ferent T and extracted the temperature dependent potential. The temperature dependent potential 
is always larger than the singlet free energy, approaching it from above. We do no see screening 
effects present in the potential for T < 190MeV. The value of the extracted temperature potential 
at large distance is very close to the value of the phenomenological potential used in Ref. [JJ. 
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